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Fractal dimensions of self-avoiding walks and Ising 
high-temperature graphs in 3D conformal bootstrap 


Abstract The fractal dimensions of polymer chains and high-temperature graphs in the Ising model 
both in three dimension are determined using the conformal bootstrap applied for the continuation 
of the 0{N) models from N = 1 (Ising model) to A = 0 (polymer). The unitarity bound below 
A = 1 of the scaling dimension for the the 0(A)-symmetric-tensor develops a kink as a function of 
the fundamental field as in the case of the energy operator dimension in the Ising model. Although 
this kink structure becomes less pronounced as A tends to zero, an emerging asymmetric minimum 
in the current central charge Cj can be used to locate the CFT. It is pointed out that certain level 
degeneracies at the 0{N) CFT should induce these singular shapes of the unitarity bounds. As an 
application to the quantum and classical spin systems, we also predict critical exponents associated 
with the Af = 1 supersymmetry, which could be relevant for locating the correspoinding fixed point in 
the phase diagram. 

Keywords Conformal field theory • Stochastic Loewner evolution • Self-avoiding walk • 0(n) model • 
3D Ising model • N=1 supersymmetry • Unitarity bound • Current central charge 


1 Introduction 

Conformal field theory (CFT) is an indispensable framework in deepening our understanding on the 
universality class of the critical phenomena which goes beyond the renormalization group (RG). Despite 
its incomparable success in 2D, the clues for 3D CFT has been scarce until recently. The recent 
breakthrough came from numerical determination for the 3D Ising exponents m using the crossing- 
symmetry sum rule for the Z 2 -symmetric intermediate states in the four-point function {4'4’4’4’) of 
same scalar field (the fundamental field ip in the A(/)'^-theory). The key empirical observation, which 
becomes a cornerstone in this so-called conformal bootstrap approach, was that the scaling dimensions 
of the spin and energy operator in the Ising model corresponds to a “kink” that emerges along the 
unitarity upper-bound curve for the dimension A ^2 of the leading non-trivial Z 2 symmetric operator 
e =: (j)^ : as a function of the dimension of the fundamental field. This singular shape, the kink in 
the unitarity bound, is shared also in the case of the sum rule for the 0(A)-symmetry, which can be 
used to map the critical 0{N) model with A = 2, 3, • • • ,00 on the A^-As plane [3] , where A^ and 
As respectively stand for the dimension of the fundamental field (pa (“a” is an 0{N) label) and the 
dimension of the energy operator e = : {(pa)'^ ■, which is the leading non-trivial operator in the 

0{N) singlet sector S. 

Towards an analytic understanding on the consequence of the 3D conformal symmetry, it would 
be important to aim at a representation theory of the spectrum generating algebra analogous to 
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the degenerate representation in the Virasoro algebra [5]. Another outstanding direction would be to 
generalize the ideas in the stochastic Loewner evolution (SLE) so as to describe critical geometry in 
3D. In this respect, the importance of the continuous family of the critical 0{N) models below N = 2 
could be emphasized more since in 2D they precisely represent the continuous family of the models 
described by the SLE^ with 2 ^ k ^ 4 via the trigonometric relation 

K = A-k! arccos(—A^/2). (1) 

Mathematicians proved that the Hausdorff dimension of SLE^ curves is given by 

dpin) = 1 + k / 8 , ( 2 ) 

(Beffara’s theorem [?])■ In physics, the same fractal dimension can be computed from the dimension 
of the 2-leg operator (a special case of the watermelon operator for an arbitrary number of legs 0 ) 
represented by an 0{N) symmetric tensor operator which behaves as a scalar under the spatial 
0(2) rotations. 

In this paper, we study the 3D 0{N) model below N = 2 with a focus on the fractal dimension of 
the loops in the high-temperature expansion. The fractal dimension may be given hy dp = D — Ap, 
where Ap is the scaling dimension 0 of the most relevant operator in the 0{N) symmetric tensor 
sector T. Apart from the models with N ^ 2, where Ap has been estimated [3], there are several 
important physical cases in 3D, where understanding based on the conformal symmetry, in particular, 
the determination of the fractal dimensions may be interesting. 

(a) Polymer {N = 0) The N ^ 0 limit of the 0{N) symmetry, where the degeneracy of Ap and 
As occurs, describes polymer chains under the excluded volume limit (a self avoiding walk) as shown 
in the celebrated work by de Gennes [5] . A direct approach to this polymer limit makes various OPE 
coefficients singular and makes the current bootstrap method, which hinges on the positivity of the 
squared OPE coefficients (a main part of the unitarity), difficult. For instance, the square of the OPE 
coefficient XP’ for the stress-energy tensor may have a simple pole at A^ = 0 since the Ward identity 
tells us that it is inverse proportional to the central charge Cp, which is essentially proportional 
to the number of components N. Practically, as N tends to zero, this pole seems to result in an 
effective slowdown of the convergence to the optimal unitarity bound, meaning that the number of 
derivatives necessary to attain a given precision increases more rapidly. Accordingly, the detection 
of the kink at N = 0.1 within a limited computational cost becomes much more difficult compared 
with the case of finite N (e.g. N = 2). We circumvent this difficulty by assuming that a clear change 
of the slope dCj/dA^ for current central charge Cj defined through the conserved current [TU] 
may correspond to the dimension A^ of the GET. This analysis leads to the estimate for the fractal 
dimensions dp = 3 — ^^(O) ~ 1.701. 

(b) The Ising model {N = 1) and the J\f = 1 SUSY fixed point The operator content of the 0{N) 
model at Af = 1 contains that of the Ising model as its singlet sub-sector and its thermodynamical 
exponents can be determined from the well-studied dimensions A^p for the spin operator and Apj 2 = As 
for the energy operator. Since the 0(1) model contains only one component scalar, it is less noticeable 
that the dimension Ap of the symmetric “tensor” ipab may carry important information on the critical 
exponents. It is, however, natural to consider that Ap is one of the geometric exponents that determines 
the fractal dimension dp = 3—A^)!) ~ 1.734 for the high-temperture graphs also measured by a Monte 
Carlo (MC) simulation for the 3D Ising model [IT]. As a natural extension of this analysis, we also 
give the fractal dimension which would possibly be realized by the magnetic flux loops in an effective 
gauge theory, which appears, for instance, in the Kitaev model plus the local exchange interaction of 
the Ising type m- An interesting possibility is that the phase diagram of this model, “magnetic three- 
state of matter” (or its slight extension), may contain the fixed point of the 3D J\f =1 superconformal 
field theory (SCFT), whose 2D counterpart is a well-established SOFT [13], which may explain the 
Majorana fermion nature of the 2D Ising model as the Nambu-Goldstone fermions associated with a 
spontaneous breaking of the supersymmetry. Unlike in 2D, where the J\f =1 SCFT corresponds to the 
universality class of the Ising tricritical point m, our view is that the SCFT and the Ising tricritical 

^ This dimension of the relevant operator in T-sector is denoted by Ap following the convention in [4]. This 
“T” should not be confused with the stress-energy tensor , which has spin 2 with 0{D) and the fixed 
scaling dimension D. 
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point are distinct fixed points in 3D. This 3D SOFT is also proposed as a boundary effective theory 
for the topological snperconductor M- 

(c) The model at TV = —2 and its possible relation to the loop-erased random walk The 0{N) model 
at iV = — 2 may be considered as an endpoint of the continuous family of the 0{N) model in the sense 
that the dimension for the fundamental field and energy operator reduce to the mean field values 0 
(Z\^, Z\s)=(l/2,1) and may be paired with the other end point N = oo, where the mean field value of 
(Z\^, At) = (1/2,1) and the spherical model value of As = 2 are realized. Among these operators in 
N = —2 and N = oo, the only nontrivial dimension is At{—2); as in 2D [57], it would be natural to 
conjecture that dp = ^ — At(—2) ~ 1.614 is the fractal dimension of the loop-erased random walk. 

Apart from the MC simluations already mentioned, there are still vast works of related simulations, 
among which some notable are certain sophisticated tests of the conformal invariance in the 3D self¬ 
avoiding walk {N —>• 0) [15], the worm algorithm that can be applied for continuous values of ^ 0 
m, and certain clever algorithms with analysis that get rid of the correction-to-scaling to attain ever 
improving precision on the self-avoiding m and loop-erased [TB] random walks. 

Our emphasis is not on the precision for the critical exponents, though some of them including 
perhaps the anomalous dimensions slightly above TV = 0 and the fractal dimension for the Ising model 
may already be more accurate than existing MC simulations WH- Instead, it is our purpose here to 
consider how the conformal invariance may be used to determine the fractal dimension, without any 
use of machine generated random numbers, and to help opening a way to understand more theoretical 
aspects (such as the kink formation, the representation theory, the 3D SLE, and so on) in the 3D 0{N) 
CFT in general. 

This paper is organized as follows. We consider the 0{N) model for a global range of —2 ^ TV ^ oo in 
Section lTm and show that the fractal dimension dp can be regarded as a geometric RG eigenvalue given 
by the dimension Ap of the traceless symmetric tensor ip at- Section [2.2l is a quantitative discussion on 
how the gap between Ap and As closes in the polymer limit TV —0 using a simple 6-loop RG analysis 
leaving the details in Appendix. In Section 13.11 the intermediate states in the four point function is 
classified into three sectors (S: singlet, T: traceless symmetric tensor. A: antisymmetric tensor) using 
the operator product expansion (OPE) (pa x (pb of the fundamental fields. The key equation in the 
0(TV) conformal bootstrap, namely, the crossing symmetry sum rule is reviewed with a brief discussion 
of the solution manifold with regard to the unitarity bound. In Section [3.21 the definitions and useful 
1/TV-expansions of the current central charge Cj as well as those of the standard central charge Cp 
are given. The implication of the unitarity and corresponding implementation of the bootstrap, though 
being standard, are given in Section iTfil We present our main results in Section[4] We give a qualitative 
description on how an effective smoothing of the kink in Ap occurs in the polymer limit TV —>■ 0 (we 
idenfify it as a severe unitarity wall, across which the continuation of the unitarity-saturating solution 
is interrupted) and discussions on how certain level-degeneracies in the 0{N) CFT would be related to 
various singular shapes (the kinks in Ap, Cp, and in particular, Cj) of the unitarity bounds in Section 
o and Section 14.21 respectively. We determine the fractal dimensions by the conformal bootstrap 
for the polymers (TV —>■ 0) in Section HTSl and for the 3D Ising high-temperature graphs (TV = I) in 
Section IT3] We compute the fractal dimension for TV = —2 in RG and conjecture that it corresponds 
to that of the loop erased random walk in Section 14.41 In Section 14.61 we estimate the set of the scaling 
dimensions {A^, A^ 2 ) for the J\f = 1 SOFT and discuss the relation to the critical exponent v as well 
as the fractal dimension dp\suSY for the corresponding excitation. We conclude with selected future 
directions in Section [SJ 


2 The 0{N) CFT for a global range of —2 ^ TV ^ C50 

2.1 The fractal dimension and the traceless symmetric tensor (pab 

We start with the discussion on the two relevant operators e and ipab in the 0{N) model, which 
respectively belong to the 0(TV)-singlet sector (S) and the 0(TV)-symmetric tensor sector (T). These 
operators are formed as a bilinear of the fundamental field (pa with the scaling dimension A^j,, which 

^ The non-renormalization of these dimensions is due to the topological property of the 0(TV)-vertex: 6cd- 

{SabScd + SacSbd + SadSbc) = (TV + 2)5ab, which vauishes at A = —2, regardless of the space dimension D |26 |. 
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transforms as a fundamental representation of the 0{N) group: 

N 

S : e{x) = (3) 

a—1 

c N 

T ■■ ^ab{x) =: (j)a(l)b ■■ -(fl- ■ (4) 

c=l 


The energy operator e is already in the single component model and plays an essential role in the initial 
formulation of the conformal bootstrap for the Ising model, which has the Z 2 = 0(1) symmetry [T]. 
The most relevant operator tpab in T sector is responsible for the crossover phenomena with respect to 
the symmetry breaking 0{N) —>■ 0{M) x 0{N — M) with an arbitrary M. As in the general statistical 
model, the two-point function {4>a{x)4)b{y)) can be expressed as a sum over self-interacting random 
walks between x and y [TH] (also The Hausdorff dimension of this random walk is given by 

dp = (i>2lv [ 23 ]: where v and cj)2 are respectively the correlation length exponent and the crossover 
exponent [22] of the 0{N) model. Since these two independent exponents are related to the scaling 
dimensions As of e and Ap of tpab by 


1 

" D-As 


4>2 = 


D — Ap 
D-As’ 


(5) 


one has a simpler expression for the Hausdorff dimension 


dp = D — Ap^ (6) 

which may be viewed as a geometric RG eigenvalue yc in the light of the fact [24] that the magnetic 
(thermal) RG eigenvalues can be determined by the relation yu = D — A^ {yp = D — Z\s). 

In the range — 2 ^ TV ^ 00 , this dimension Ap decreases monotonically from a certain value 
Ap{—2) (FIG. m and below for the meaning) to the trivial value Ap{qo) = D — 2 and at TV = 0 
crosses the dimension As of the energy operator e (in the 0{N) singlet sector S, as mentioned), 
which in turn increases monotonically from As{—2) = D — 2 to As{oo) = 2 in the same range of 
TV (one may also notice the asymptotic slopes computed in 1/TV-expansions are symmetric in 3D: 
-g^As{oo) = —-Q^Ap{oo) = 8). Actually, this somewhat dual behavior of Ap and As in the global 
range of TV is almost shared in the 2D 0{N) model though the range TV G [—2, 00 ] should be replaced 
by TV G [—2, 2], where the model has a critical point and exact results are available from the Goulomb 
gas m, SLE [B], and GFT torus partition funtion (as described just below) for continuous values of 
TV; it is also likely to be a generic feature of the 0{N) GFT in 2 ^ D < 4 from the RG point of view. 

The operator content of the 2D 0{N) model with —2 ^ TV ^ 2 can be studied exactly by the torus 
partition function [29] . Using the Goulomb gas coupling g {1 ^ g ^ 2) determined by the relation 
TV = —2cos(7r(;), the central charge (which is a 2D counterpart of Cp in Section |32|) and the scaling 
dimensions are given by 

^^ = l_g/2_3/(8g), As=4/g-2, Z\t = 1 - I/(2g). (7) 

Since the SLE parameter k is actually related to g by k = 4/g, the last relation in (0 with (|6|) is 
equivalent to the formula 0 in the Beffara’s theorem |7|. 

In the 2D torus partition function, the multiplicity TV(TV -|- l)/2 — I=(TV — 1)(TV -|- 2)/2 for the 
traceless symmetric tensor ipab tends to zero as TV —>■ 1 in accordance with the observation that the 
expression in (|4|) apparently vanishes at TV = 1. It is, however, instructive to note that the dimension 
Ap = 5/8 in the TV = 1 model {g = 4/3) is of physical relevance. Namely, it corresponds via (l6|) to 
the fractal dimension dp = 11/8 of the Ising interfaces, which are the SLEk =3 curves. The relevance 
of this tensor ipab for generic D ^ 2 in the 0{N) sum rule at TV = 1 will be discussed in the end of 
Section ED and will be used to determine its scaling dimension in D = 3 in Section 14.51 Similarly for 
the 2D TV = —2 model, the dimension Ap = 3/4 {g = 2) leads to the fractal dimension dp = 5/4 of 
the loop-erased random walks (SLEk =2 curves) |2Z] • 

We will give in Section a simple estimate for dp in the TV = — 2 model using m by a pseudo 
e expansion in the 6-loop RG, which agrees with the numerical simulations results obtained for the 
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Fig. 1 The scaling dimensions of the singlet scalar e (As: solid red) and traceless symmetric tensor (fiat (At: 
dashed blue) as a function of in 2D (left: eq.ljT])) and in 3D (right). The right branch of the unitariy bound 
(dashed gray) with the Z 2 case 32] is shown for 2D as a guide to the eye. The 3D curve is obtained as the 
[5/l]-Pade approximant for N € (—2, 7] continued by the curve from the pseudo e-series (Appendix) for N > 7 
and should be regarded as schematic as the anomalous dimension tends to be smaller than the genuine value. 
The N ^ 00 asymptotics (dotted) are shown for both scaling dimensions. 


3D loop-erased random walk We use the conformal bootstrap to determine At in the 

0(1) model and the fractal dimension dp of the high-temperture graphs in the 3D Ising model [TT] in 
Section 14.51 


2.2 The degeneracy of the relevant operators from S and T sectors in the limit N ^ 0 


In addition to the above two important cases, we are especially interested in the iV —>■ 0 limit of the 
3D 0(N) model, which describes dilute solutions of polymers, where the random walk becomes self¬ 
avoiding. Besides such physical relevance, the limit TV —>■ 0 is theoretically special for the following two 
reasons. First, some of squared OPE coefficients may become negative for TV < 0 due to single poles 
at TV = 0, which makes it difhcult to take the approaches based on the unitarity, on which most of the 
present conformal bootstrap schemes depend. Second, as mentined above, TV = 0 is the precisely the 
point where the degeneracy of the two scaling dimensions As and At take place. As a quick example 
using d?]) in 2D, As = At = 2/3 (dp = 4/3) 0 follows from g = 3/2 for TV = 0 and the gap opens with 
the following asymmetric TV-derivatives, 


d 

ON 


As = 8/(97r) = 0.282942-•• 

7V=0 


d 

dN 


At = -l/(97r) = -0.035367 • 

N=0 


( 8 ) 


It would be notable that this derivative for At for 2D is somehow almost unchanged in magnitude for 
3D as we will see below. The leading term in the e-expansion may be compared with ([5|) as 


d(As — At) 


dN 


N=0 


I/tt D = 2, 

e/8 + 0(e^) D = 4-e, 


(9) 


where each contribution of the derivative for D = 4 — e is d(As, AT)/dN = (3e/32, —e/32). In Ap¬ 
pendix, we compute a pseudo e series using the input of the 6-loops D = 3 RG calculations m and 
present reasonable estimates by a simple Fade analysis together with the best-known results of the 
e-expansion up to e®. As a simple estimate, we take the average of the six and five-loops and the 
maximum deviation as an error. This gives, 


d 

dN 


As 

N=0 


0.1238(28), 


d 

dN 


At 

N=0 


-0.036(7). 


( 10 ) 


® It is also well known in 2D that dp = 4/3 coincides with the Flory value [2^ dp = v ^ (D + 2)/3. 
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The same analysis for the derivatives at = 1 (the Ising point) yields 


d 

dN 


As 

N=1 


0.1017(35), 


d 

dN 


Nj' 

N=1 


-0.032(5), 


( 11 ) 


The derivatives for At^ which are useful in this paper, increases only slightly (~ 10%) in the interval 
N G [0,1]. Nevertheless, in order to get better estimates, one may found it more useful to keep both 
(na and m than to choose one of these two. More concretely, the variation At{N 2 ) — At{Ni) with 
0 ^ iVi < iV 2 ^ 1 can be better approximated by {N 2 — Ni) times the derivative at the midpoint 
(iV2+iVi)/2 obtained as a linear interpolation between (TTOl) and dill) . For instance, a roughest estimate 
for the variation between iV = 1 and TV = 0 may be obtained as /1 t(1) — At{0) = (+1) x (—0.032(5) — 
0.036(7))/2 = —0.034(4), where the errors in (fTUl) and (fTTl) are assumed to be independent. Although 
we do not use the last example, which would maximize the uncertainty, one may check 0 this estimate 
may reasonably connect the results obtained independently by conformal bootstrap in Section 14.31 
(iV = 0) and in Section IT3] (N = 1). 


3 Operator product expansion of the fundamental fields in the 0{N) CFT 


3.1 Crossing symmetry sum rule 


The crossing symmetry sum rule used in this paper is the most basic one (in the sense it does not involve 
the mixed correlators [3T]) in the conformal bootstrap for the CFT with a global 0{N) symmetry [H 
131131133 as described briefly below. The fundamental field in this theory is a scalar operator (pa, 
which transforms as an 0(iV)-vector, with dimension A^. Crucially, the OPE of (pa with itself may be 
decomposed into three sectors; 


pa{x) X (pb{0) 




^ab 


E 

(A,e)&T 


AA,£'-^(ab) + 


Z 

(A,e)GA 


^A,£^[ab] ’ 


( 12 ) 


where S, T, and A denote the 0{N) singlets sector of even spin, the 0{N) symmetric tensor sector of 
even spin, and the 0{N) anti-symmetric tensor sector of odd spin, respectively. Note that there is an 
infinite tower of the scaling dimensions {A^^\ A^p, • • • } for the states with fixed f in each sector. The 
dependence on x is omitted on the right hand side. The sets of the OPE coefficients 
(with X = S, T, A and the tensor labels are omitted) encode important dynamical information in the 
0{N) CFT and satisfy highly nontrivial constraints due to the associativity of the operator algebra. 
These constraints can be expressed as a sum rule that follows from the equivalence (the crossing 
symmetry) of the two different expansions of a single four point function {pa{xi)pb{x 2 ) 4 >c{x 3 )(pd{x 4 )) 
from the two distinct degeneration limits (xi -A X 2 and xi -A X 4 ), where the contribution from the 
identity operator, which belongs to the singlet (S) sector, becomes dominant. 

Let us write the contribution from each sector in the OPE (TT^ in the channel xi -A X 2 as follows: 


S' = 


>x,eGAAu,v), >^A,eGA,iA(A, A = 

{A,t)eS' {A,£)€T 


^A,e^A,£{u,v). 

{A,e)GA 


(13) 


Here in the first sum, the set S' includes all the operator in ^-sector except the the identity operator 
{A,£) = (0,0), for which the contribution for the four-point function is simply {xl 2 X 34 )~‘^‘'’AbPcd, 
which is usualy the dominant contribution in the limit xi -A- X 2 - For concreteness, we note that the 
conformal partial wave (global conformal blocks) GA,i(u,v) in H-dimensions is a function of the two 
cross-ratios given as 


u = zz = 


2 2 
^12^34 
2 2 ' 
^13^24 


V = {1 — z)(l — z) = 


rjA 

X14X23 
2 2 ’ 
2 ^ 13^24 


X'i'j — Xi X i 


(14) 


This paper focuses on the approach to A = 0 through At- However, some cru des t benchmark for the 
singlet dimension is possible using As(l) = 1.41264(6) obtained for the Ising model [30] as follows: A s(0) ~ 
1.41264(6) + (—1) X (0.1238(28) -I-0.1017(35))/2 = 1.300(22), which is consistent with At( 0) in Section 14.31 as 
expected. 









7 


and has the following form [36] in terms of the radial coordinates re*® = 2/(1 + 

00 

GAAr,0) = ^5^BSr^+"C';(cos0), (15) 

n=0 j 


/’£'\ l-l-f_ 

where the coefficients i?„ t with j = i + n,i + n — 2, - ■ ■ ,iiiax{i — n, ——^-) can be iteratively fixed 

by the Casimir differential equation for the H-dimensional conformal group, and CJ with i' = (£> —2)/2 
is the Gegenbauer polynomial. 

The crossing symmetry can be nicely seen using the OPE (fT^ and the notation (TT^ as follows: 


(t>a{xi) MxA 
4>b{x2) 4 >c{x3) 


{xl2xlA‘^* 




{)( (l + 5') + 

(1 + 5') 


+x-^)( 


r- 


A 


)(+X- 


N 


T- 


)(-X 


A 


(16) 

(17) 


where all the possible three tensor structures after the contractions are represented as )(= SabAd, 
X= 5ac5bd, SadSbc, and the tilde notation is used for representing the quantities with u and v 
interchanged, thus S' = S'{u,v), S' = S'{v,u), etc.. Using the notation X± = v~^*X ± u~^*X 
{X = 1,5', T, A), the 0{N) sum rule follows by comparing the terms for each tensor structure )(, 
and X: 


5(_ + ^1 - —^ 71 + .4_ — -1_, (18) 

5+ - ^1 + —^ T+ - A+ = -1+, (19) 

r--A-= 0 . ( 20 ) 

For each iV, the solution manifold for the crossing symmetry (HU)-® consisting of the points repre¬ 
sented by the effective GET data (the possible set of scaling dimensions and spins {A,£) in X-sector 
with associated OPE coefficients A^^) may be inifinite dimensional. An important one-parameter- 
family solution, which is conveniently parametrized by A^, can be singled out along the boundary 
of the unitarity (dictated by the lower bounds (l30l) and the positivity (|31l) '). whose projection onto 
A^-At plane is shown for each N in FIG. |2| As is well-known, the search for this unitarity saturating 
solution can be formulated as a linear optimization problem (see Section [3.3l for more details) and can 
be solved with the aid of knowledge on the global conformal blocks Gzi_^(it,?;) such as that in (fTSl) . 

In order to clear up a common source of confusion, it is worth to make a careful distinction between 
the spectrum of the 0(1) model {N = 1) and that of the Ising model in our formulation. The Ising 
(Z 2 ) sum rule, which is used in [1] for instance, follows from the crossing symmetry of the four-point 
function for a single scalar {<f){xi)4>(x2)4>{x?,)4‘{xA))■ Since only the singlet fields appear in the Ising 
OPE 4>y~ (j), the whole contribution in xi —>■ X 2 except that from the identity operator may be denoted 
as 5'. Then the crossing symmetry leads to 1 + 5' = {u/vA*{l + 5'), which further simplifies to, 

S'_=-l-. (21) 

Now the logic is as follows. A proper subset (ITS)) and (EUl) of the 0{N) sum rule for X = 1 implies 
the Ising sum rule (l2T|) . which, along with the requirements of saturating unitarity (' (l30ll and OT), 
is sufficient for a given A^ to single out a unique solution for 5'. Thus, in particular, the unitarity- 
saturating solution for the Ising sum rule may be embedded into the solution for the 0(1) sum rule as 
its 5-sector. In this case, one may still generalize this Ising spectrum to a solution for the 0(1) sum 
rule, which may also admit non-empty T 0 A sectors, in addition to 5-sector, determined in turn by 
solving —37+ — A+ = —1+ — 5(^ and 71 — A- = 0, where the Ising contribution 5(|_ may be regarded as 
a seed generating these sectors. In particular, T-sector in this solution contains the rank-2 symmetric 
tensor operator ipab in (|1|) with a non-vanishing squared OPE coefficient also for X = 1, which actually 
determines the fractal dimensions for the Ising high temperature graphs as shown in Section 14.51 















3.2 The central charges Ct and Cj 

If we omit RG irrelevant operators and keep only most important ones, the OPE (TT^ becomes 


</'o X '5^6 — (l + ^As,0 


D,2 


T^'") Sab + J[ab] - 


( 22 ) 


where and J^j, are the stress-energy tensor and the conserved vector current, respectively. As 
the symbol in (1^51) signifies, is a spin-2 0{N) singlet with dimension D and is a spin-1 
anti-symmetric vector (J^j, = — with dimension D — 1, which transform as an 0(iV)-adjoint 
representation. The conformally invariant two-point functions of T^^{x) and Jabi^) m, 

^ Jr,(x2)) = [SaeSbd - SaJbc), (23) 

®12 

with normalization given by the surface of unit (D — l)-sphere Sd = jr{D 12) and with 

/M^,P^(^) = 1 (^IPP{x)I'''^{x) + I^'^{x)I''P{x)^ - IP''{x) = SP'' - 2^^. (24) 

The Ward identities for the stress-energy tensor T and the conserved current J leads to 


NCT/CT-,lree = ' , Cj/Cj;free = 1/ (A^_i,i) ' , 


(25) 


where CT,free = ND/{D — 1) and Cyfree = 2./{D — 2) are free field values. For the 0{N) GET, some 
useful results are known in the IR fixed point. These include the e = A — D expansion 


r ir = 1 - 5 N + 2 

t ! T;free 12 (iV-f 8)^ 

and the 1 /A^-expansion in I? = 3 0 

CT/CT-ireeO{N) = ^ ^ ' 


0{e% 


o 


1 


r 'c - 1 3 fv -I- 2 2 


0(e3), 


CjICj-hee = 1 — 


64 1 
Qtt^ N 


o 


1 


(26) 


(27) 


Regarding (IMl) and (EZl), it is an open direction to study nnunii under which conditions these Ct and 
Cj in 3D are monotonically decreasing along the RG as the central charge does so in the 2D unitary 
system [35]. Together with the leading correction in 1/iV in the expansion [42] . 

1 _ 256 32 (3^^(-3402C(3) - 61 + 1081og(2)) - 3188) 1 

2^ 21tt^N^ 2437r6iV3 ^ N^’ ^ ^ 

the large N asymptotics of the central charges Ct and Cj may be written as a function of as 

10 S2 

CT/C'T;free = 1 " -J^ i^l) ’ C'j = 2 " y + O {^1) , (29) 


where the slope dCr jdA^ has been already considered in |3] . On the other hand, we observe that the 
other one dCj jdA^ shows interesting behavior (FIG. |3|) for the case iV ^ 0 <C 1 of our main interest 
0. The formation of an effective minimum near iV ^ 0 is discussed in Section 14.21 and used to estimate 
the fractal dimension dp oi polymer chains in Section [4.31 


® For 1/N coefficient, there is a mismatch by a factor 2 between (4.25) and (6.8) of [TOj. Our results on the 
slope (I29II in FIG. [3] as well as mm supports the value —64/97r^ in (I27II reproduced from (4.25). Also by 
using a Fade analysis on (I28II . the universal curves A^-As and A^-At in FIG. [T] can be drawn, which will be 
discussed elsewhere. 

® We thank Tomoki Ohts uki for pointing out that similar kinks in Cj can be observed via the direct Cj 
minimization in D = 3 | 43| for N ^ 2, and Yu Nakayama for further discussions. For Ct in 3D, it is known 
that the direct Ct minimization reproduces Ct along the unitarity bound (via As-maximization) for A = 1, 
but not for generic A > 2 [4]. It is possible to check similar characteristics are shared by Cj in 3D. It would be 
also interesting to study the implication of these phenomena on the solution space of the crossing symmetry. 
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3.3 Conformal bootstrap, the unitarity bound, and the search space 


A particularly important one-parameter-family solution of the crossing symmetry (Ull)-® lies along 
the boundary of unitarity, which may connect the whole spectrum of the free theory and the 0{N) 
CFT. This solution can be singled out by a linear optimization as mentioned in the end of Section [3.II 
The unitarity consists of the following two conditions. First, the scaling dimensions in a D-dimensional 
theory must satisfy the lower bounds, which correspond to the requiment that the anomalous dimen¬ 
sions be positive) [HIH^HHIITT] . 


A ^ 



+ D-2 


for 

for 


0 , 

£ = 0 , 


(30) 


where the inequalities are saturated by conserved currents such as and {£ > 0) and by free 
scalars such as a fundamental field (pa at the free theory (£ = 0). Second, the squared OPE coefficients 
must be positive: 

operators in (IT^ . (31) 

These two requirements of the unitarity enable one to solve the crossing symmetry di)-® along the 
unitarity bound via the simplex algorithm or the semi-definite program m- 

We use the standard simplex algorithm (Sec. 6 of [5]) with our particular implementation based 
on the code [48) . As usual, the simplex algorithm is used in order to try to determine if there exists 
a solution of the crossing symmetry (HD)-® that satisfies the lower bounds (1301) and the positivity 
for a given and in the region Ax > Ax,, for the dimension Axg of some low-lying operator; 
here we use it for the symmetric tensor (pab {Ax = At) mainly for the reason given in Section [4.11 6. If 
solutions do not exist (do exist), the next search region At > At^ can be chosen narrower such that 
At-i^ < Ato {At-^ > Atq). Then one may take a bisection procedure from some initial finite interval of 
At and narrow the search region by each trial. If the initial interval is taken wide enough, the iteration 
eventually reaches the upper-bound At^ for At, for which the solution in At ^ Ay^ is expected 
to be unique (in particular. At = At^). The value of At^ is measured numerically by setting the 
bisection accuracy goal 6{At), which we typically take S{At) less than 10“"^. 

In practice, the crossing symmetry constraints are extracted by a truncated Taylor-expansion 
around the symmetric point u = v = 1/4 of the sum rule (I18p - (l20p also with a truncated number 
(■■max of the spin sectors. The simplex algorithm (at j-th step of the bisection) searches the spectrum 
region bounded from below by (l30ll (with At ^ Ajv.i) and from above by an appropriate upper bound 
Amax, which should be taken large enough. The derivatives of (TTRll - d^nil are computed with respect to 
the coordinate (a, 6) defined from {z,z) = (a -I- \/b,a — \/b)l2, which is related to the cross-ratios as 
(u, v) = (zz, (1 — z){l — z)). Following the convention in [35], they are reduced to the following set of 
the derivatives of the global conformal block GA,i{a-, b) in (flKll . which one may select as 


|9™9^Gzi,f(a = 1,6 = 0) m = 0, • • • , 2(ninax - n)-b m^ax; n = 0, • • • , nmaxj (32) 


with some (mmax, ?T-max), which consists of/C = (minax + ?T-max +l)('ramax +1) derivatives. In general, the 
unitarity bound becomes more strict for larger number K, of the derivatives, although an exact form of 
convergence to the optimal bound is not well-understood so far. Also for a given /Cq, the bound usually 
depends only scarcely to a paticular choice of (mmax, ?T-max) with K, ~ /Cq- The number of spins .^max 
should be taken large enough with respect to the choice (mmax, ''T-max) so that resulting bound does not 
depend on ^max- Our default choice used for measuring the critical exponents is (mmax, »^max) = (8,8) 
with /C = 153, Amax = 70, and (max ~ 50. 


4 Singular shapes of the unitarity bounds 

4.1 Smoothing of the kink in the unitarity bound for the tensor dimension At 

A crucial observation that gave an initial momentum to the recent revival of the bootstrap studies was 
perhaps the emergence of the singular shape (a kink) at A^ = A^; ising in the unitarity bound curve 
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Fig. 2 The unitarity upper bound for the dimension At of the tensor ipab obtained with the derivative cut 
oft set at (nmaxiUimax) = (8,8) for N = 0.1 and (umax, TUmux) = (7,2) otherwise. We have also verified that a 
higher cut off makes the kinks sharper for generic N, but not so effectively for N ^ 0.2. 


of the scaling dimension A ^2 = A^ 2 {Atj,) in the case of Z 2 symmetry studied for the Ising model Q . 
Actually, the appearance of singular shapes in the conformal bootstrap seems more ubiquitous as we 
will see also in this work. 

Below we use the sum rule (TTSll - dSnil for non-integer N. It is important to realize that the unitarity 
is violated in the free 0{N) models for any non-integer N [M], and it is also most likely to be so in the 
IR fixed point of the 3D 0{N) model for non-integer N. In particular, one should not expect the IR 
fixed point lies exactly at the unitarity saturating solutions. Nevertheless, the results seem to point to a 
weak effect of violation that the IR fixed point may still lie close to the unitarity-saturating solution as 
long as N is positive and is not too close to N = Q. For TV < 0, a severer effect of the violation prevents 
one from obtaining meaningful bounds as discussed in the end of Section 14.1.11 More discussions on 
the effect of the unitarity violation for generic non-integer TV > 0 is postponed to Section [5] 

We here show the unitarity upper-bound obtained by the bisection for At (Section 13.311 in the 
0{N) model with 0.1 ^ iV ^ 2 in FIG. [2] This result may be regarded as an extention to the region 
iV < 2 of the bounds previously obtained for At ior N ^ 2 (and As for ^ 1) [4], where observed 
single kink has been used to estimate the scaling dimensions {A^,Ay) in the 0{N) CFT for each N 
and for each sector Y = S,T. As one can see, the curves in FIG. [2]for iV ^ 0.5 has a clearer kink and 
these can be used to estimate {A^,At) of the 0{N) GFT; note that once A^ is determined, other 
low-lying dimensions are also available since the simplex algorithm has reached a unique solution (see 
the end of Section [XT]) for the crossing symmetry (IT^ - (I^ along the unitarity bound for a given A^. 

Remarkably, however, the kink becomes less and less pronounced as N tends to zero (the polymer 
limit) thus practically making the determination of A^ more difficult. Even in such circumstances, 
other singular shapes may remain in other universal quantities along the unitarity bound. This is 
indeed the case, and we observed a very clear change of the slope dCj/dA^ in the current central 
charge Cj as discussed in Section 14.21 and used this for the limit —>■ 0 in Section 14.31 Now two 
important remarks are in order. 


^ In view of the picture that there are infinitely many Z 2 symmetric primary operators above e =\ : whose 

levels are separated by non-trivial intervals and may be repulsive to each other, the observed straightness of 
the lowest level A ^2 on the right side of the kink {A^ > A^-^ ising) is also remarkable. 
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4-. 1.1 Convergence to the optimal bound. 

The first is on the convergence of the bound to the optimal shape with respect to the truncation (1321) 
of the derivative orders (m,n) on the conformal block 6), which is currently unavoidable in 

numerics. Although the upper-bound obtained in a finite truncation is rigorous, a larger number K, of 
derivatives leads to more restrictive bound (i.e. a lower upper-bound) and actually makes the kink more 
sharp, where the convergence to the optimal bound tends to be faster than that in the rest. In addition 
to this, the convergence becomes much slower in the polymer limit iV —i 0. Thus a bruteforce approach 
to the limit is taking K, large enough with respect to a given small N. In practice, the unitarity bound 
for At looks smooth for N < 0.2, which makes the detection of the kink (a discontinuity in the slope) 
becomes very hard within a reasonablely large number of derivatives (/C = 153). Here it would be also 
worth noting the present guess on the optimal shape (/C = oo) for the limit —>■ 0. On the right of the 
kink, the convergence of the slope dCrjdA^ to a almost constant, as in the cuves for iV ^ 0.5, seems 
to be plausible. On the left, finite-/C curves are convex upward as in FIG. [2l It seems, however, there 
are no strong indications that excludes the possiblity that the optimal shape on the left is also almost 
straight in 3D, while there is a 2D example where the bound is likely to be convex downward |49j . In 
general, it would be interesting to consider if there is a principle that forbids the optimal unitiarity 
bound for these scaling dimensions be convex upward with some reasonable assumptions. We will give 
a qualitative argument in Section 14.21 on the enhanced slope BAt / dA^ on the left of the kink in view 
of the level dynamics. 

The slowdown of the convergence is probably related to the degeneracy of two levels At = Ag and 
the severe unitarity wall at = 0, where some squared OPE coefficients 0 including 2 ^ for the 
stress-energy tensor has a pole, which can be seen by the relation (1^51) from the Ward identity and by 
our observation that the ratio Ct jCT-Svee remains finite ~ 0.955. The detailed analysis on this special 
limit > 0 and quantitative knowledge on the order of convergence in the conformal bootstrap in 
general might be useful and would deserve further investigation. 

4 . 1.2 Gap assumptions. 

The second remark is especially relevant if one tries to find the upper bound for As by the Ag-bisection 
in the case A < I. In [3], it was found useful to complement the unitarity lower bounds (1301) by an 
extra gap assumption that scalar fields {i = 0) in the right hand side of the OPE (fT^ be bounded 
below by the RG canonical dimension: 

A ^ D — 2 for£ = 0 operators in the R.H.S. of the OPE (fT^ . (33) 

This has previously been used in order just to improve numerical stability [4]; in particular, the 
condition At ^ 1 in D = 3 was not supposed to change the resulting solution of (HE])-® through the 
Ag-bisection. Indeed, it can be checked that this gap assumption makes no distinction in the resulting 
spectrum for A ^ 2. 

In the bootstrap for A < I, however, we found that the Ag-bisection using the pure unitarity 
condition (1501) and (1511) may yield a solution that violates this additional gap assumption (1551) , namely, 
a solution with 1/2 < At < 1 < Ag, which satisfies (1501) . but can not correspond to the 0(A) GET 
from the RG point of view. This phenomenon may remind us of a level repulsion between At and Ag 
in the solution space that becomes stronger as A —>■ 0. Accordingly, if one sticks to keep the extra 
condition At ^ I, the Ag-bisection yields another unphysical solution 0 which contains At = 1 < Ag 
and thus saturates the extra gap assumption set by hand. In contrast to the Ag-bisection, our bisection 
for At, which should be the lowest dimension scalar contained in the product 4>a x (/& in the 0{N) GET 

® Another important example of a singlular OPE coeffici ent in A —^ 0 may be for three energy operators, 
which would play a role in the mixed correlator bootstrap m- The physical origin of the divergenc e of Xlg can 
be traced back to the strong repulsion between the loop segments e in the 0{N 0) loop model | 50| . 

® Around A = 1 the saturation of dSSl) may not be so serious as the kink in As appears around the 
expected Ising position, which should be consistent with our observation that the OPE coefficient XJq of the 
(unphysical) level At = 1 is negligible compared to those for other operators. However, below A = 1 this makes 
much difference: for instance, a kink in Ag emerges even in A = 0.1, which was smoothed in the solution with 
the pure unitarity conditions. Again, it is obvious that this solution with At = 1 can not represent a physical 
spectrum. 
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with > 0, is free from such problems. In particular, the extra assumption Ag ^ 1 does not change 
the resulting solution. As we have just seen, besides the direct role of determining the fractal dimension 
(p]). the property that At has the lowest dimensions (like a ground state in quantum mechanics) in 
the right hand side of the OPE (fT^ adds the study of At a special importance. 


4.2 The slope change in Cj and the level degeneracy in A-sector 

The singular shape is not restricted to the unitarity bound for scaling dimensions; it may also appear 
in the unitariy bound for the OPE coefficients of the conserved currents such as and which 
are via the Ward identities reflected on sudden changes of the slopes in the central charge Ct and 
the current central charge Cj defined in (l23|) . We show CT/{NCfree) and Cj for 0.1 < < 2 in 

EIG. [21 most of which has just one point where the slope change occurs. As the interaction in the 
0{N) GET becomes infinitesimal in the N ^ oo limit, both changes in {A^,Ct) and {A^,Cj) from 
the free field values tends to zero. This corresponds to the asymptotics of these central charges given 
by (IMl); for a finite N, not too small, these slopes are actually shared as the initial slopes for small 
A^ — 1/2, where the effective interaction may be weak. In the case of our interest {N —>■ 0), a change 
of the slope dCjjdA^ is enhanced and an effective minimum is formed for N < 0.5. The effective 
minimum is used to estimate the dimension A^ of the fundamental field in A^ —>■ 0 in Section 14.31 
Before this application, we present some preliminary analysis on the mechanism (a level dynamics 
along the unitarity saturating solution that connects the free field theory and the 0{N) GET, where 
the levels are essentially eigenvalues of the infinite dimensional matrix obtained by linearizing the RG 
flow around the fixed point in the theory space) behind the formation of these kinks. 

As A^ = 1/2 corresponds to the free field theory, the effective anomalous dimension rj = 2{A^ — l/2) 
may be considered as an effective interaction parameter. If one traces the spectrum of scaling dimensions 
along the unitarity bound, one will meet a reorganization of the spectrum when crosses the value at 
the kink A^*, which is expected to be A^ of the 0{N) GET as a similar phenomenon has been observed 
in the Ising model [2]. Along the unitarity-saturating solution for the crossing symmetry of the 0{N) 
sum rule (UHl)-®, this reorganization may be qualitatively different for N > 0.5 and N <C 0.5, which 
may be described as follows. Suppose we superpose the curves for the sub-leading scaling dimensions 
on FIG. m and a certain dimension bifurcates to the right (left) as we move to larger (smaller) 
then let us call this i?-bifurcation (L-bifurcation). 

Before describing more details of the 0{N) spectrum, let us mention that the usage of the word 
“bifurcation” here does not necessarily mean the bifurcation of the common square-root type. We take 
the following important case to illustrate this. On the Ising spectrum obtained via the Z 2 sum-rule 
(EU, an extensive description is given in Section 3 of Ref. [2], where the recombination is shown to 
occur a la Hilbert’s infinite hotel. As both sides of the IR fixed point (A^ = A^*) have infinitely 
many operators, the effective correspondence between the two spectrum can be non-trivial depending 
on the versions of the infinite hotels, for which “00 = 00 -I- 1” and “00 = 2 • 00 ” are described below. 
In particular, if one temporarily denotes the scalar operators {£ = 0) on the right (A^ > A^*) by 
and those on the left (A^ < A^*) by £,x,£',£" ■ ■ ■ in ascending order of the scaling 
dimensions, the recombination of the spectrum is observed as in the following: 

■■■ , £'" A£'', £” A£', £'Ax, £^£- (“oo = oo-bI”), (34) 

where the symbol “A’’’ stands for a connection of the nearest levels by a sudden descent from the right 
of A^ = A^* to the left, and x is a decoupling (null) operator, which only appears numerically on the 
left side (A^ < A^*) with a small (ideally vanishing) squread OPE coefficient. Although the nearest 
levels in this case (e.g. £” and £') never touch with each other within a finite numerical bootstrap, it 
is also suggested that the recombination transition becomes shaper and eventually the nearest levels 
would be connected in the limit of the infinite number of derivatives (/C —> 00 ). In this regard, the 
connection between both sides for a large, but finite /C would also seem to be, 

... ^ A £"' -b £" {^£' -b x), A £" -b £' (-bx), £' A £' + x, f (“00 = 2 • 00 ”), (35) 

where the levels on the right bifurcate (or possibly multifurcate) to the left as represented by the 
symbol Here in the connection £" A £" -b £', for instance, the lower branch £" —>■ £' causes 
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a sudden decent by a finite gap > 0, while the upper branch £" —^ £" (eventually 

connected as —>■ —0 in /C —?► c») does not meet such a large jump if /C is large enough. 

It is also worth to mention that if the decoupling operator x on the left indeed disappears in the ideal 
limit /C —> oo, the connection is then more like 

• • • , £'" £"', £" £\ £' ^£', £^ £, (36) 

where all the connections of the levels may be continuous (with no gap left) and may have smoother 
(or even no) kinks compared with that of the lowest one {£). In the scalar (^ = 0 in S-) sector along 
the unitarity-saturating solutions for the 0{N) sum rule (TTSll - d^ni) around TV = 1, we observe a similar 
recombination as in (1341) as expected. Besides more common looking bifurcation, below let us simply 

call these processes in (1551) e.g. £” £" + a L-bifurcation of £". 

For TV > 0.5, we observe that an i?-bifurcation of the dimension of the sub-leading spin- 1 

antisymmetric tensor in ^-sector 0 (just above the conserved current J^j,) and a L-bifurcations of the 
sub-leading dimension = A^i ~ 3.8 from S'-sector[3 occur simultaneously at A^ = A^* of the 
0{N) CFT (see (1151) for the sectors S, T, and A). After the L-bifurcation, the lower branch of ^ 
flows into the free value (A^, A) = (1/2, 2) (being RG unstable, the free theory at A^ = 1/2 tends to 

have fast varying subleading dimensions; thus it is numerically subtle to see lim^^_,,i /2 ~ 

seem to contribute to larger slopes dAs/dA^ and dAx/dA^ in A^ < A^* via the level repulsion. 

Now for N ^ 0.5, the i?-bifurcation does no longer coincide with the L-bifurcation; the latter may 
be observed at much larger A^. The lower branch after the L-bifurcation of the subleading spin-1 
dimension flows into the level Zl = 2 of the conserved current J^f, just below in the same A-sector. This 
isolated L-bifurcation and the confluent behavior in A-sector should lead to the enhanced change of 
the slope dCj/dA^ (i.e. the effective asymmetric minimum) at Zl^ = Zi^* of the 0{N) CFT and the 
subsequent divergent behavior of Cj in Zl^ > Zl^*, respectively. The change in the slope dCrldA^ 
seems to be mainly due to the L-bifurcation in the L-sector; the L-bifurcation in A-sector may also 
change dCx jdA^^ but this effect is rather weak as in the curve for TV = 0.1 in FIG. [S] The discussion 
above is obviously not enough to fully describe the level dynamics as one moves along the unitarity 
bound and to understand how it may lead to the formation of the kink. On the other hand, the 
degeneracy of the levels at A^ = Zl^* might play a role in constructing the putative representation 
theory for 3D CFTs. Therefore, the bifurcations and the confluent behavior observed here may deserve 
further studies. 


4.3 Determination of A^ and the fractal dimension in the limit TV —>■ 0 

There are practically at least three ways to estimate A^ (or equivalently, the anomalous dimension r\) 
of the 0{N) CFT with some reasonable assumptions for each case: 

1. Calculate A^ that gives the asymmetric minimum of the current central charge Cj (FIG. [3]), 

2. Calculate A^ where the L-bifurcation of the subleading spin-I dimension occurs, 

3. Locate the kink {A^,At) in the unitarity upperbound of At (FIG.[2). 

As discussed in the previous section, the method 1 and method 2 are essentially equivalent and should 
give consistent estimates. Within our derivative truncation (1551) of /C ~ 153, the method 1 is applicable 
for TV < 0.4. Although the method 3, which will be used for TV = I in Section 15751 has an advantage of 
giving simultaneous estimates for {A^,At), it may not be so accurate for TV < 0.2 as the smoothing 

An anal ogou s R-bifurcation of a spin-1 operator is also observed in the Af = 2 supersymmetric (SUSY) 
Ising model [52] , where the decoupling operator of the lower branch never touches the level of the J/), at A = 2. 
It is also remarkable that the N = 0 model has a tw isted TV = 2 SUSY in 2D [55], whose origin, the presence 
of underlying Osp{2M, 2M) for any TU in A = 0 |54| . is actually independent of the space dimension D. 

This L-bifurcation of A^^ (dimension for £' =: : with E = 0^) is accompanied by a level crossing 

of Ay ^ (dimension for : EFab ■ with Fab = (j>a4’b — E/N) and Ay ^ for A ^ 1. In the Ising model, Ag^^ becomes 
A^4 ~ 3.8, which gives the correction to scaling exponent cu ~ 0.8. In the XY model (A = 2), we reproduce 

Ay^ ~ 3.65 [57][5n][55]. More detailed study of the subleading spectrum is beyond the scope of this work. 
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Fig. 3 The current central charge Cj (left) and the central charge Ct (right) as functions of A^. The curves 
are obtained with the derivative order cut off at (rimax, m-max) = (8, 8) for N — 0.1 and at (umax, rUmax) = (7, 2) 
otherwise. The large N asymptotics (I2S1) with the slopes —32/3 and —10/3 are also shown (dashed). 


of the kink inevitably occurs (Section 14.II) . Taking these into account, using the criterions 1 and 2 for 
a zoom-up data, one may obtain by the conformal bootstrap 

= 0.5145(2) for N = 0.1, (37) 

which amounts, by an ad-hoc linear extrapolation from = 0.518151(6) for the Ising model {N = 1) 
[50] . to = 0.5141±0.0002 for = 0 with uncertainty simply copied from (1571) as we will not use 

this later (the notation ±a;, instead of (x), is temporally used to indicate uncertainties for the ease of 
comparison). There are various RG estimates for the anomalous dimension ry = 2Zl0-|-2 —H (see |60] for 
review), which have in general relatively larger uncertainties than those for the exponent = D—Ag. 
Nicely, the estimate for = 0 above lies almost at the center of the result = 0.5142 ± 0.0013 

obtained from updated and the most accurate RG computation [^. Note also that MG simulations 
are rarely able to measure this with an exception of = 0.5125 ± 0.0007 m. which is slightly 

smaller than our bootstrap result. 

Now, we estimate the optimal unitarity bound for At at N = 0.1 with in (1371) and extrapolate 
it to A^ = 0. The unitarity upper bound A^ for At approaches the optimal value from above as the 
number of derivatives K, tends to infinity as shown in Table [TJ The downward uncertainty for the last 
digit (10“®) of each Atp{lC) due to the choice of the bisection accuracy is shown as a subscript. Also 
note that these last digits may be subject to change by the choice of the cut-off for spins. As already 
mentioned, the convergence for this small value of N becomes much slower than that for generic N like 
N = 1. Note also that At*{IC) has a minor variation in addition to the overall tendency to decrease 
(the upper-bound must decrease as the constraints gets stronger). This is expected since At* depends 
on the precise choice of the derivatives (I32L which has more information than just a one number /C. 
The uncertainty induced by this variation, however, does not become dominant in the analysis below. 

The estimate for the optimal bound At*{oo) in Table [1] is obtained by a phenomenological fit 
using [3 At*{JC) = At*{oo) + -|- 0( ^pVi ) p = 2, which is presumably better than the 

raw bound AT*{K.^ax) with /Cmax = 153 obtained by the bisection. We adopt this value At*{oo) = 
1.2948(36) as optimal for N = 0.1 with this conservative error bar, which includes the entire residual 
\At*{JC ms,x) — ziT*(oo)| ^ 0.0029 as well as the smaller uncertainty propagated from (1571) along the 
curve in FIG. [2j Using the above value At*{qo) for At at N = 0.1 and the derivatives BAt/ON in 
([TU|)-(fTT|) in the extrapolation from N = 0.1, one may have. 

At = 1.2984(36) for A^ = 0, (38) 

where the error due to the extrapolation is estimated as the uncertainty in m multiplied by 0.1 giving 
7 X 10“^, which is negligible compared to the other uncertainties. We note that the extrapolations from 
other small values of N would give consistent estimates meaning that the RG extrapolation by OAt /dN 


It is equivalent to find the intercept at a; = 0 in the linear fit for the data {x,y) = (1//C^, At*(/C)). 
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IC 

80 

99 

117 

135 

153 {ICmaA 

OO 

A^ 

At*{N = 0.1) 
At* (TV = 1) 

1.30517s 

1.266272 

1.300964 

1.265982 

1.300384 

1.265942 

1.297994 

1.265903 

1.297705 

1.265893 

1.2948 

1.2654 

0.5145 

0.51815 


Table 1 The unitarity bound for At with N = 0.1 at = 0.5145 and with = 1 at A^ = 0.51815 ob¬ 
tained with various choices (rimax, mmax)=(7, 2 ) 30 , ( 8 , 2 ) 99 , (8, 4 ) 117 , ( 8 , 6 ) 135 , ( 8 , 8)153 of the derivatives, where 
the number K, of derivatives is give n as a subscript. The data is not ideally smooth as it depends weakly on 
the precise choice of derivatives (EH), which is not completely specified by 1C. For fixed IC, the bound At* may 
be smaller by up to y x 10 “®, where y is shown as the subscript of the last digit. 


is correct and actually avoidable in principle. With the degeneracy As = At at iV = 0 understood 
f Section 12.2L this value (1551) is consistent with As = 1.2999(32) (y = 0.5882(11)) from the RG [5T| 
and with As = 1.29815(2) [v = 0.587597(11)) from the most accurate MC, which is far ahead of other 
simulations in accuracy Using the relation ([5]), the symmetric tensor dimension (1551) leads to the 
fractal dimension. 


= 1.7016(36). (39) 

Here let us just mention that this is much larger than the Flory value dp = = 5/3 = 1.6666 • • • 

[25) and should be more precise. The comparison with the corresponding results for As{= At) from 
more modern litertures are just given below (I38L 

Now, two remarks are in order. First, we note that using another choice p = 1 in the fit would 
lead to a larger residual |Z\T*(/Cmax) — 4iT*(oo)| ~ 0.0085 resulting in At = 1.293(9) for N = 0. 
Although this value from p = 1 is still consistent with the other estimates, one clearly needs to further 
increase /Cmax from /Cmax = 153 in order to obtain a better estimate, which is computationally time 
consuming. Although there seems to be no decisive difference between p = 1 and p = 2 regarding 
the quality of the hts, the ht using the function with the coexisting powers p = 1 and p = 2 yields 
\At* iJCraa.^) — 2iT*(oo)| ^ 0.0024 giving At = 1.2988(32) for N = 0, which is effectively the same 
result as (|55|) obtained with p = 2 only. Second, we comment on the subtlety in the analysis for 
non-integer N. Before doing so, let us briefly summarise on the three quantities all for N = 0.1: 
4iT*(Aniax), At*{oo), At, for which the level of rigor is decreasing in this order. The unitarity bound 
2lT*(Aniax) is a I'igorous upper bound, albeit not optimal. The extrapolation to the optimal bound 
At* {oo) involves the phenomenologica fit, for which the entire residual \At* (ICma-x) — 4\t*(oo)| is 
included as an error; the resulting estimate is not rigorous, but would be called conservative. Last but 
not the least, we use At*{oo) as an estimate for At in the 0{N) CFT. This is justified if the unitarity 
bound is saturated at the 0{N) CFT, which is emplically expected at TV € N as will be done for TV = 1 
in Section 1431 On the other hand, one should not expect an exact saturation for TV = 0.1 as suggested 
by the unitarity violation in the free 0{N) model for non-integer TV [M]. We nevertheless expect that 
\At — At*[oo)\/At 1, which means that the unitarity-saturating solution still passes very close 
to the location of the 0{N) CFT as it happens in the bootstrap for the Ising model in non-integer 
dimensions H ^ N [3] , which is shown to be non-unitary |81j . This point is corroborated in Section [5] 
Our purpose here is not a pursuit on the numerical precision, which is presently less than the 
MC [17) . but is giving a new perspective on how the conformal bootstrap can be used to determine 
the fractal dimension. Although more sophisticated approaches to the limit TV = 0 deserves further 
consideration as in Section jS) the result here may be already encouraging enough to let us believe 
gaining a deeper understanding on the 3D self-avoiding walk based on the conformal invariance is 
promising. 


4.4 The end point TV = —2 and the loop erased random walk 

The present form of the conformal bootstrap, which depends on the unitarity, can not directly be 
applied to the 0(—2) model, despite its importance as an endpoint of the continuous 0{N) family that 
would be paired with the spherical model limit TV = oo ('Section l2.1ll . Instead, we compute the pseudo 
e expansion (r-series) for the symmetric tensor dimension At in the 3D 0(—2) model from the result 
of the fixed dimension 6-loop RG |5T1[55] . More details can be found in Appendix, where our parallel 
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1.66667 

1.59259 

1.62462 

1.60200 

*1.64174 

*1.57244 

1 

1.71429 

1.57143 

1.61495 

1.61136 

1.61642 

1.61649 

- 

2 

1.62406 

1.62277 

1.61106 

1.61385 

1.61649jjQg] 

- 

- 

3 

1.62279 

1.62418 

1.61643 

1.61622 

- 

- 

- 

4 

1.60818 

1.61698 

1.6162[i4] 

- 

- 

- 

- 

5 

• 1.63063 p,] 

1.61633 

- 

- 

- 

- 

- 

6 

*1.59218 

- 

- 

- 

“ 

„ 

- 


Table 2 The Fade table for dF|jv =_2 = ^ — At- The positive real pole closest to 1 is shown in the bracket. 


analysis for the iV-derivatives of At and As is performed. The result is, 

3-Z\t= 2-^ - 0.0740741r2 + 0.0320229r3 - 0.0226127 t^ + 0.0397418t® - 0.0693066r® + O ( t ^) 

O 


The 5- and 6-loop simple Fade analysis parallel to that in Appendix using Table [2] yields, 

dF\N^-2 = 1-614(16) (41) 

with which the relatively large uncertainty comes from the oscillating data along the boundary of Table 
[2] (i.e. the direct series M = 0 and its dual L = 0) and is estimated as a root-mean-square deviation. 
As expected, the estimate dni) is consistent with the fractal dimensions in the loop-erased random 
walk (LERW) obtained in a number of works including c^lerw = 1.614(11) by the functional RG [62] 
and by various simulations: c?lerw = 1.623(11) [23]) 1.6183(4) [23], 1.6236(4) [25], 1.62400(5) [TSj. As 
a remark on the analysis of the Fade table, it may be possible to note that field theories tend to prefer 
a slightly smaller central values compared with the numerical predictions. In our case, omissions of 
the four boundry data (indicated by * in Table|2|) would lead to dF|jv =-2 = 1.6162(23), which would 
agree with the functional RG m and older simulations PlIM]. but would be dehnitely smaller than 
the most recent numerical results |231[I1]. It would be interesting to improve the present conformal 
bootstrap so as to analyse N < 0 across the severe unitarity wall, to obtain a better estimate on 
dFl 7 v=- 2 ) dimension for sub-leading operators that is responsible for the correction to scaling. 

Such study may contribute to deeper understanding on the LERW from the conformal invariance in 
3D. 


4.5 Fractal dimension of the high-temperature graphs in the 3D Ising model 

The fractal dimension of the critical excitation in the 0(1) model is most straightforwardly and rig¬ 
orously accessible by the conformal bootstrap since the bisection for At yields a clear kink along the 
unitarity bound just as in the conventional cases W ^ 2 |4| and since the 0(1) model is unitary so 
that the 0(1) IR fixed point may be expected to saturate the optimal unitarity bound as in the Ising 
case. By a brief inspection of the zoom-up of FIG. H] one obtains At ^ 1.266 and A^ ^ 0.5181, for 
which the latter is consistent with the known result for the Ising spin operator A^ = 0.518151(6) [30] . 
It is also empirically interesting, though being not rigorous at all, that if we adopt an ad-hoc criterion 
that the kink is located around where the bisection takes longer than a certain period (^ 2 weeks, 
for instance), we obtain another precise estimate = 0.518149(6) from the 0{N) sum rule (fT31) - (l3(Il) 
with = 1; near the kink, the optimization takes indeed longer time than the other generic points, 
while the convergence to the optimal bound with respect to the number of derivatives K, becomes much 
faster. 

The unitarity upper bound At* for A = 1 at A^ = 0.51815 for various choices of the derivatives, 
represented by its number of the components /C, are shown in Tabled) This leads via ([21) to an estimate 
for the fractal dimension for the high-temperature graphs in the 3D Ising model, 

dFU^i = 1.7346(5), (42) 

which agrees well with dp = 1.7349(65) from the 3D Ising plaquette-update MC simulation [TT] and 
with the worm algorithm simulation dp = 1.734(4) [T2]. Our results by the conformal bootstrap (A = 0 










Fig. 4 The scaling dimensions As (solid) and At = D — dp (dashed) as a function of N in 2D (left: exact) 
and in 3D (right: [4/l]-Pade). Only for N = —2 (LERW) in 3D, the error bar is visible where At is obtained by 
the RG. For = 0 (SAW) and = 1 (Ising high-T graph), At are obtained by the 3D conformal bootstrap 
where the error bars are less than the line width. 


(l39)l and A^ = 1 (l42)l ) and by the RG {N = —2 (IdlTl ) are shown in FIG. [4l The related estimates using 
(ini) for the XY model {N = 2) is given in a footnote in Section!^ 

Since the 0(1) model is unitary, the rigorous upper-bound Z\T*(/Cniax) applies to the 0(1) GFT, 
which yields dp = 3 — Z\T*(Afmax) > 1.7341. In contrast to the case with N = 0.1 in Section H31 
it would be natural to assume that the 0(1) GFT saturates the optimal bound. Note also that the 
uncertainty 6 x 10“® in as being multiplied by {dAT/dA^)\A^=o. 5 i 8 i 5 , may induce the error only 
less than 10“^ in dp. Then the source of uncertainty is virtually restricted to the extrapolation of 
At* from /C = /Cmax to /C = oo. This uncertainty |zlT*(Afmax) — 21 t*(oo)| ^ 5 x 10“"^ as shown in 
is evaluated from the fit using At*{K.) = At*{oo) + + O(^c^) with p = \. We keep this 

conservative estimate, while the uncertainty is roughly halved if one uses the fit with p = 2 or becomes 
even smaller if one uses an exponential fit, for which the quality of the latter seems decent for this 
particular case. It would be very useful if one had a general theory for the scaling of the residual with 
respect to the optimal bound as a function of /C. In any case, the conformal bootstrap here for this 
particular fractal dimension seems to give an order of magnitude more precise result compared with 
the MG simulations mM- 


4.6 The tricrtical Ising fixed point and the M =1 SOFT 

In 2D, it is well-known that the GFT for the trictitical Ising model (c = 7/10) is known to be the first 
member of the series of the minimal J\f =1 superconformal CFTs (SCFTs) [13]. In condensed matter, 
the N =1 SOFT in 3D is propsed to describe boundary excitations in the topological superconductor 
[14] . The action of the one-component Gross-Neveu-Yukawa model in D-dimensions is given by 

S = j d^x (43) 

where the Yukawa coupling mixes the bosons (j) and the fermion tl). At the Af = 1 supersymmet¬ 
ric IR fixed point, the interaction is described in terms of the superpotential W = and a real- 
supermultiplet E = (j) + Oij) + 0"^ (jp [66], where 0 is a fermionic coordinate (in superspace) of dimension 
— 1/2 regardless of the space dimension D and 9^ = e“^0Q0/3 (a, /3 = 1,2). Thus it is pointed out that 
the intersection of the following extra constraint 

A^r^=A^-^ + l (44) 

with the unitarity bound curve A^i = A^ 2 {A(j,) may be used to locate the SOFT [5S] . 
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1.3537(2) 

0.5705(2) 
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1.3541(2) 

0.5713(2) 

3.0443(27) 

1.3542(2) 

0.5722(4) 

3.0323(51) 

1.3532(4) 

0.5731(6) 

3.0235(80) 

1.3520(6) 

0.5737(5) 

2.9955(64) 

1.3528(5) 


Table 3 The cut-off JC of derivatives and the corresponding position of the intersection between the 

SUSY relation ^ and the linear fit ((451) for the unitarity bound for A ^2 as well as the coefficients a and b. 


In contrast to the branch of this curve in the left of the kink, the right branch in the relevant region, 
which intersects with da, is observed to be almost linear both in 2D and 3D [5]. Let us then define 
the coefficients a and b of the linear fit in 3D; 




.2 



(45) 


For each choice of the set of the derivatives (n-max, m-max) = (2fc, 1), (fc = 4, 5, • • • , 10) with the cut-off 
K. for the number of derivatives given in Section [331 and the spin cut-off at L = 50 for k = 9,10 and at 
L = 40 for otherwise, the linear fit is performed for the right branch of the the unitarity upper-bound of 
A ^2 in the range G [0.54, 0.59] obtained through the Z 2 sum rule. TableOshows the position of the 
intersection as well as the coefficients of the fit El with respect to 1C. Although the convergence 

of the data may be improved by choosing the cut-offs (and the range of the fit) which consumes much 
more computation time, this would not seriously change the qualitative argument below. Table [3] may 
lead to an estimate 


0.574, + 1, (46) 

which is slightly larger than the one-loop RG result 1/2-1-1/14 = 0.571 • • • [2] and satisfies the rough 
lower bound Z\^=i > 0.565 [66] B ■ 

The geometric exponent for the Af = 1 fixed point may be obtained by the conformal bootstrap 
giving the symemtric tensor dimension At 1.43 in the 0(1) model at A^ = The derivation of 

dni) does not seem to depend on whether a fixed point may allow the description by the supersymmetry 
or not, although more detailed analysis would be useful. By assuming ([51) holds here also, one would 
obtain. 


^f\susy ^ 1-57. (47) 

In order to find the fractal object with the dimension (14711 in a lattice model, one natural candidate 
would be the Kitaev model augmented by a local Ising exchange interaction at hnite temperature, in 
which the non-supersymmetric fractal dimension has already been realized by the magnetic flux loops 
m- In particular, the critical temperature along the paramagnetic to quantum-spin-liquid transition 
was determined by identifying the magnetic flux loop in this extended Kitaev model (an effective Z 2 
gauge system) with the 3D Ising high-temperature graphs and by using the knowledge [TT] of the 
fractal dimension for the latter, for which our independent estimate is in (1421) . Similarly, it would be 
interesting if the fractal dimension (SB is realized by the flux loop at some point in the phase diagram 
(or its suitable extension) and can be used to locate the emergent Af = 1 fixed point by simulation. 
Indeed, the phase diagram is already very rich and studied for understanding the effect of thermal 
agitations on the topological order, and more specifically, the thermal fractionalization of the quantum 
spins into Majorana fermions. It is also worthwhile to note that if the Af = 1 SOFT is realized in the 
vicinity of the “tricritical” point B in [l2]; the corresponding RG fixed point may be different from that 

It would be interesting to study this slope a as a function of D < 4, though it is beyond the scope here (see 
[3] for related figures). One ad-hoc interpolation between the two values a ~ 3 (D = 3) and o = 8/3 {D — 2) 
is o = 2(6 — D)/{5 — D) inspired by the approximate coincidence of a with the critical dimensions where (ff 
and 0® becomes RG marginal. 

The same approximate value A^f, = 0.565 has recently been reproduced by the fermion bootstrap m- 

Numerical accuracy in the Monte-Carlo simulation m has still to be improved to measure the exponents 
at this partic ular tricritical point. We thank the author, Y. Kato for useful discussions and clarification on the 
related work mi- 
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of the 3D Ising tricritical point, which is widely believed to be described by the mean field exponents 
in view of the RG argument including (jfi interaction |()9j . 

Another interesting model which would be relevant to the Af = 1 supersymmetry is the Blume- 
Capel model [72] • In 2D, the tricritical point of this model belongs to the tricritical Ising universality 
class, which can be identified with the celebrated Af = 1 fixed point [13] ■ In the latter, the set of 
scaling dimensions which is obtained as an intersection between the relation (1441) and the unitarity 
bound (which is expected to be saturated by the analytic solution [ 68 j to the crossing symmetry along 
A ^2 = |z \0 + |) is actually = (1/5, 6/5), where the energy operator e and subleading 

energy operator e' both in the Z 2 -even sector together form one operator in the Neveu-Schwarz sector 
(while the spin operator, which is Z 2 -odd, independently belongs to the Ramond sector with dimension 
Aa = 3/40 [I3])- Thus the correspondence (A^, A^ 2 ) —>• {A^, Ag) in the Ising model should be replaced 
by A-^ 2 ~^) —t (4\e, A^r) in the Af = 1 model in 2D. Regarding the crossing symmetry, this can 

happen in 2D because cr in the Ising model and e in the Af = 1 model have the same fusion rule since 
they belong to the same position (r, s) = (1,2) in the Kac table of the Virasoro representation. If we 
use the same identification A^^^ —>■ Ag as in 2D, this gives the exponents in the thermal sector 

2-a ~ 1.236, i/-0.412, (48) 

which are in a reasonable agreement with the exponent 2 — a = 1.213 found by the variational RG 
method in the phase diagram of the 3D Blume-Gapel model m- As noticed in [73], the latter value, 
which was aimed at the 3D tricritical exponent, deviates considerably from the standard Ising tricritical 
exponent, which is believed to be the mean field value 2 — a = 3/2 |69j as already mentioned, while 
their results agree impressively with the 2D exact results for both the tricritical and the critical point, 
and with the modern estimate for the 3D critical point. This seems to leave some possibility that their 
method actually detects an additional Af = 1 fixed point with the exponent (H51) in the Blume-Gapel 
model other than the possibility that their value quoted above is simply a very poor estimate for the 
3D Ising tricritical point. In any case, in view of the RG, it is likely that it is only in D = 2 that the 
Ising tricritical point and the Af = 1 fixed point can be identified with each other. It is interesting to 
study how these two fixed points would deviate from each other in D = 2 + e [zniEi] and evolve all 
the way to the different universality classes in D = 3 by the conformal bootstrap. 


5 Conclusion 

We take the simplest bootstrap approach based only on the crossing symmetry of the four-point 
function {(t)a4'b<t>c4>d) of the fundamental fields for the one-parameter-family of the 3D 0{N) model 
with a special focus on the fractal dimension dp for the range 0 ^ A^ ^ 1. Besides the property of 
being exactly at the severe unitarity wall (Section |4T]), the limit A^ —>■ 0 may be characterized by the 
degeneracy of the two operator dimensions As and Ap in any space dimensionality D. Accordingly, 
a more elaborate approach to such limits would need to deal with the possible logarithms that could 
appear in the four-point functions [THlITTirTH] . Also if one tries to perform the mixed correlator bootstrap 
[3T1I32] including the energy fields e, for which the logarithms appear at the level of two point functions, 
it would be inevitable to face with these logarithms. In that case, the smooth continuation by the 
conformal bootstrap to the 2D problem from D = 2 + e would be also interesting since the correlation 
function including e in the 2D 0{N) model can be dealt with both in the integral-representation 
[7^ and the differential equation using the degenerate representation [5] at the integer level 3 in the 
Virasoro algebra. This is a subject of further research. 

The issue of unitarity violation is certainly important and nontrivial, and it is currently not obvious 
to what extent it affects the estimates in Section l473l for non-integer N obtained using the assumption 
of the unitarity, which is actually weakly violated as outlined below. Nevertheless, there seems to 
be various possible improvements for the study of the bootstrap for non-unitary systems (e.g. the 

Another identification ^ A^ yields v ~ 0.701, which seem to agree with four estimates for v 

(~ 0.71) by the functional RG [74] for the Af = 1 UV Lagrangian. Note that the latter u is not meant for 
a physical realization (e.g. the Blume-Capel model) and is just an indication that computed from the 

mass renormalization may agree with the conformal bootstrap. In the 2D Af — 1 fixed point, for instance, the 
observed value 2 — a = dv = 10/9 ~ 1.11 (73] follows from A^^^ = 1/5, but not from A’^'T^. 
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determinant method [82] or the extremal function flow method [83]) suggested from the present work. 
To elaborate on some of these, let us first recall that the positivity is violated in the free 0{N) 
models at any non-integer value of N, which can be shown by computing the norms □ for a certain 
class of composite operators [51]. Although the above argument only applies to the free model, the 
recent study on the unitarity violation in the single scalar ())'^-theory {N = 1 is fixed) in the fractional 
dimensions D = i — e points to a quatitatively similar behavior, both for the free (UV) theory and 
for the interacting (IR) theory, that the positivity violation indeed occurs, but generically only at the 
operator with very high scaling dimensions D [5T] , which would explain why the bootstrap for the 
Ising model in non-integer dimensions works decently [3]. Our exact computation for the squared OPE 
coefficients in the 2D 0{N) model for non-integer N also shows a similar behavior for the positivity 
violation. In view of these analysis, we conclude that the set of the scaling dimensions {A^,At) of 
the IR 0{N) FP for non-integer N should not saturate the unitarity upper bound At ^ A^{A^) in 
general, but hopefully it lies in the vicinity of the bound (it could appear on either side of the bound), 
for N not too close to = 0, at which the squared OPE coefficient for the low-lying operator, namely, 
the stress-energy tensor {A = D) changes its sign due to the pole as discussed in Section 14.1.11 
For N ^ 0, on the other hand, if the family of the unitarity saturating solution obtained here along 
At = A^{A^) passes indeed nearby the IR 0{N) FP, an inclusion of one more spin-1 operator just 
above the conserved current in the truncated spectrum in the determinant method |82] should 
qualitatively improve the behavior of the non-unitary solution according to the mechanism in Section 
14.21 for the emergence of the kink in the current central charge Cj. It would be nice to quantify the 
effect of the unitarity violation in more detail and clarify how (via the flow method |83] for instance) 
the putative non-unitary IR 0{N) FP and the unitary solution here could be connected with each 
other. 

In the long run, it would be interesting to generalize the ideas in the SLE so as to describe the 
critical geometry embedded in 3D though the task would severely face to our limited understanding 
on the 3D geometry since the SLE is intrinsically based on the Riemann mapping (uniformization) 
theorem on the 2D conformal map. In this respect, the 3D 0{N) model with —2 ^ TV ^ oo offers a 
natural one-parameter-family of the loop ensembles m for an extensive simulations), which is likely 
to have a conformally invariant measure m for a simulation at TV = 0). We have been able to focus 
on the region 0 ^ TV ^1, which have the most interesting cases EE as the two boundary points (namely, 
the self-avoiding-walks in the TV —>■ 0 model and the high-temperature graphs in 3D Ising model in 
the = 1 model) and to estimate the fractal dimension dp = ^ — At by the conformal bootstrap 
using the unitarity conditions. Although the positivity (13111 prevents us from applying the conformal 
bootstrap across the severe unitarity wall at A^ = 0, we compute the fractal dimension for A^ = — 2 
by the 6-loop RG giving dp ^ 1.614, which would be encouraging to conjecture that the 3D N = —2 
model may also describe the loop erased random walk as in 2D although some elaborate operator 
correspondence may be necessary in view of the logarithmic corrections m- It is of interest to see if 
some generalization of the Beffara’s theorem ([2]) exists in 3D and if the loop ensemble in the 3D 0{N) 
model with —2 ^ A^ ^ oo of the fractal dimension 1.614 < dp ^ 2 (parallel to 5/4 ^ ^ 3/2 in 

2D) has a natural parametrization in terms of the inverse-trigonometric ([T]), or other transcendental 
function of N, as k of the 2D SLE^. 
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These norms are considered in the AdS 4 /CFT 3 holography for large N. We thank Tassos Petkou for 
pointing out that there would be an interesting physics also for small N, where the bulk d egre es of freedom 
in this higher-spin holography may be related with a polymer-like product of the singletons [SS] in the replica 
limit. 

The same idea applies to the case A > 1. For some A ^ 2, At has been determined from the conformal 
bootstrap [4|. For instance, we will have dp = 1.76437(108) f or t he high-tem pera ture graphs in the XY model 
(A = 2), which agrees with the simulations dp = 1.7626(66) [TT], 1.7655(20) [5D], and 1.765(3) |16| . 
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Table 4 The Fade table for the derivative dAs{0)/dN. The positive real pole closest to 1 is shown in the 
bracket. 
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Appendix A Computations on the Af-derivatives by the fixed dimension RG 

We use the fixed dimension RG {D = 3) augmented by the pseudo-e series [86] to circumvent the accumulation 
of the intermediate systematic error due to the determination of the coupling g* at the non-trivial fixed point. 
The beta function is generalized by a new parameter r such that at r = 1 it reduces to the original beta 
function f5{g) in 3D: 


P{9,r) = -Tg + I32{g), 


(49) 


where /? 2 (p) = Pig) + g starts at order g^ with a positive coefficient of order 1. Then the critical exponents can 
be expanded in r lay eliminating g by using g = (;*(r) which solves P{g* ,t) = 0. We compute the r-series for 
the derivatives of As and At at the degeneration point [N = 0) based on the six-loop 3D RG results for Pig), 
V, 7"' m and r/T = g + 4’2/1' — 2 [55] . The results are, □ 


d 

dN 


d 

dN 


As =|rT -b 0.0241609r^ -b 0.0036762r® -b 0.0021794r‘‘ -b 0.0004817t® - 0.0019009r® + O 

= ^e + 0.0361328e^ - 0.0198967e® -b 0.03806686"^ - 0.0611648e® + O (e®) , (50) 

At = - -^r- 0.0099826t^ -b 0.0049121r® - 0.0029128r'‘ -b 0.0066431r® - 0.0101052r® + O (r’’) 

= - - 0.0146484e^ + 0.032921e® + 0.028875e'‘ -b O (e®) , (51) 


where we also show the derivatives computed from e = 4 — d expansion up to known orders |87II88| just for 
comparison. Since both r-series do not show strong asymptotic behaviors with factorial growth of coefficients 
up to the orders presented, even the naive direct summation of the series would be of some use; in particular, 
it is clearly better than the direct sum of the e-expansion. 

A simple Fade analysis, however, may improve the stability of analysis as usual. This can be irrustrated 
as follows. We show the values from the Fade approximants [M/L] for these derivatives in Table [4] and Table 
[5] respectively. The positive real poles closest to r = 1 are shown in brackets for the six-loops (anti-diagonals 


In the process, this computation naturally reproduces the r-series for 1/2 = 3 — At in [^ . 
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- 

- 
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-0.0619539p .,j. 

- 

- 

- 

“ 

- 


Table 5 The Fade table for dAT(0)/dN. The positive real pole closest to 1 is shown in the bracket. 


L + M = 6) and five-loops {L + M = 5) order approximants. For each derivative, the data occuring with a pole 
in [0.5,1.5] (indicated by *) is omitted since it is rather close to r = 1, where the series is to be evaluated. As 
a simple estimate, we tak e th e average of the six and five-loops and the maximum deviation as an error. This 
gives the value quoted in (nni and (1111) in the text. 








